The photon absorption edge in superconductors and gapped ID systems 
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Opening of a gap in the low-energy excitations spectrum affects the power-law singularity in the 
photon absorption spectrum A{yi). In the normal state, the singularity, A{Q.) oc [D/{^ — fith)]'^, 
is characterized by an interaction-dependent exponent a. On the contrary, in the supeconducting 
state the divergence, A{Q.) cx {D//S)°'{Q. — f2th)~^^'^, is interaction-independent, while threshold is 
shifted, ilth = fith + A; the "normal-metal" form of A{Q.) resumes at (f2 — Slth) ^ Aexp(l/Q). If 
the core hole is magnetic, it creates in-gap states; these states transform drastically the absorption 
edge. In addition, processes of scattering off the magnetic core hole involving spin-flip give rise to 
inelastic absorption with one or several real excited pairs in the final state, yielding a structure of 
peaks in A{Q.) at multiples of 2A above the threshold frequency. The above conclusions apply to 
a broad class of systems, e.g., Mott insulators, where a gap opens at the Fermi level due to the 
interactions. 

PACS numbers; 74.25.Gz,74.50.-l-r,73.40.Gk 
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I. INTRODUCTION 

It was demonstrated more than 40 years agoi^^i^ that 
electron a;-ray absorption coefficient in metal, A{bj), is 
strongly modified by attraction to the localized hole left 
behind. The threshold behavior of absorption coefficient 
was found to be 



Aiuo) = A ( - 



(1) 



In Eq. ^ and thereafter, uj = Q, — Jlth stands for the 
difference between the photon energy and the core-hole 
energy measured from the Fermi level, and D is the band- 
width. Prefactor, A, contains the square of the dipole 
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FIG. 1: (Color online) Schematic illustration of elastic absorp- 
tion, (a), and inelastic absorption, (b). Blue lines illustrate 
creation and annihilation of a virtual pair that participates 
in elastic absorption. Final state of inelastic absorption is 
electron with energy e and a real pair, Brown lines 

in (b): since inelastic absorption is possible only for a spin- 
ful core hole, in-gap states created by this holel4 can also 
participate in absorption. 




FIG. 2: (Color online) Absorption spectrum near the thresh- 
old for spinless (green) and spinful (red) core hole. 



matrix element between the level and the conduction 
band. In the simplest case of a weak short-range at- 
traction, V{y) < 0, of electron to the hole the expression 
for the exponent a ^ 1 has a form 



a — 2iyQ 



drV{r) 



(2) 



where vq is the density of states at the Fermi level (we 
neglect the correction, — originating from the An- 
derson orthogonality catastrophe,^). Since the diverging 
absorption Eq. ^ comes from all energy scales between 
u; and D, it is quite robust. In a finite system, the thresh- 
old behavior depends on additional energy scale, the level 
spacing**. 

Interest to the singular behavior of A{uj) near the 
threshold got a boost after it was predicted^ that this be- 
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havior manifests itself in the resonant-tunneling current- 
voltage characteristics. This prediction was later con- 
firmed in numerous exp er iment s^'^'^'^ i i o 1 1 1 1 1 ^ . Enhance- 
ment of absorption Eq. ([1]) was derived under the as- 
sumption that the density of states, v{uj), is constant 
v{lo) = within the entire frequency interval, {—D^D). 
If there is a gap, 2A, at the Fermi level the threshold be- 
havior of A(uj) is singular even without interaction with 
a hole: 



illustrated in Fig. [Jd. More additional pairs in the final 
state give rise to anomalies at uj — tOn — (2n-|-l)A, which 
have the form 



SAjuj) 
A{nA) 



2n I 1/2 nl \ 



(7) 



II. DERIVATION OF EQ. (gD 



A{lo) oc v{lo) = 



(W2 _ A2)l/2 



(3) 



and diverges near the edge of the gap. For small a it 
could be expectedi^ that this strong bare singularity is 
weakly affected by the excitonic effectsi. Indeed, the 
low-energy, < 2A, many-body processes across the gap, 
responsible for Mahan singularity, are suppressed. This 
reasoning suggests the form of the absorption in super- 
conductor 



(4) 



Eq. Q crosses over to the conventional behavior Eq. ([T]) 
at high frequencies, lo, such that Q;ln(Lj/A) ~ 1; in this 
frequency domain the effect of superconductivity is neg- 
ligible, since w ^ A. 

Even stronger modification of the absorption spectrum 
takes place, when the core hole possesses a spin, so that 
the interaction with excited electron includes exchange. 
In this case two new physical mechanisms come into play. 
Firstly, a core hole creates in-gap states^^ with binding 
energy Eq ^ a^A measured from the edges. These states, 
in turn, affect dramatically the elastic scattering of ex- 
cited electron transforming the near-gap absorption into 



A(u:) 



V2 



,1/2 



- A) + £0 



(5) 



see Fig. [2] The absorption is zero at the threshold and 
resumes (w — IS)~^I'^ falloff only for (cj — A) Eq. As a 
"compensation" of the suppressed absorption, a (5-peak 
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(6) 



emerges at the position of the bound state. 

There is another many-body feature in A(lo), which is 
specific for the exchange interaction with core hole. This 
feature originates from the fact that exchange interaction 
of electron with localized magnetic impurity in metal can 
be accompanied by creation of an electron-hole pair—. 
The underlying reason is that localized spin emerges as a 
result of the on-site Hubbard repulsion of two electrons. 
On the other hand, with electron-electron interaction, 
two electrons can be excited by a single photor^'^. In 
the presence of a rigid superconducting gap, this process 
starts from the threshold^* lo = uji = 3A, which cor- 
responds to inelastic absorption with electron and addi- 
tional pair in the final state. This process is schematically 



A. Time dependent superconducting Green 
functions 

An efficient way^ to derive Eq. ([T]) is to consider 
scattering of excited electron by a transient potential, 
V[T)9{t), and perform calculation in the time represen- 
tation. In this representation the Green function of 
the normal metal Go(t) = / dwe*"* ^/{'^ " ± «0) 
(+ or — depending on sgn(^q)) has the form Go(i) — 

—uq (t — iD~^sgn{t)) , where D is the bandwidth. 
Generalization of the scattering approach to supercon- 
ductor requires the time representation of the supercon- 
ducting single-particle Green function 



GiLO,q) 



A+(g) 



iO 



A-(g) 
+ ey-iO' 



(8) 



where = 



A2 is the spectrum of superconduc- 



tor; = VpQ with q — [k — kp) being the momentum 
measured from the Fermi momentum, k^, and Vp is the 
Fermi velocity. The projection operators ^±{q) are 2x2 
matrices 

_U 

A 2 "I" /el I A 2 \ 

(9) 
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1± 



A 



1 ± 



with following properties: Aj_(g) — A±(g) and ^.+ {q) + 
A_(g) — 1. In the basis of eigenfunctions of the 
Bogoliubov-de Gennes Hamiltonian, interaction with the 
short-range potential is described by the diagonal matrix 



V 



1 

-1 



(10) 



Note that time-dependent 2x2 Green function of a su- 
perconductor, obtained as a result of integration dwe*"* 
of Eq. ([8]), and subsequent summation over momentum, 
q, can be conveniently expressed in terms of zeroth and 
first-order Bessel functions, namely 



G{t) 



( Git) Fit) 
- [ Fit) Git) 



(11) 



where the normal and anomalous Green functions, G(i) 
and Fit), are given by 



Git) 
Fit) 



ttAj/q 



Dt>l 



Dt>l 



sgn(i) zJi(A|t|)+yi(A|t|) , (12) 



iJoiA\t\)+YoiA\t\ 



(13) 
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FIG. 3: (Color online) Conventional arrangement of timea^ 
in n- fold integral Eq. (|16p describing contribution to the 
response function due to n successive scatterings by the core 
hole. Time intervals, \ti—ti+i\, are distributed unevenly; cen- 
tral interval corresponding to line 1 is only slightly smaller 
than Remaining intervals contained in the boundary el- 
lipses are ^ Inset; blowup of the right end of the line 1. 



In the limit A — > the normal-metal Green function, 
Ga{t), is recovered from Eq. by using the small-f 

asymptote Yi{At) « -2/(7rAt), while F{t) 0. 



B. Shape of the absorption edge 

In superconductor, we generalize the response function 
to a 2 X 2 matrix, L{t), so that the absorption coefficient 
is given by the diagonal matrix element 



Aico) 



Aq 



Rc 



dt exp{—iLut) L{t) 
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(14) 



As a result of matrix generalization, the expansion of the 
response function in powers of a, 



(15) 



has the 2x2 coefficients, Ln{t), which are given by the 
following n-fold integrals^'^ of the single-particle Green 
function, G{t), 

/o M 
dh-- dt,,G{-h)VG{ti-h)V ■ -VGiU-t). 

(16) 

In the normal metal, evaluation oi A{u)) is based on exact 
analytical result^ for the infinite sum 



E 



ra=0 



21^0 



dti 



dtnGoiT-ti)---Go{tn~T') 



Go(r- 



(t-r)r' 



{t - t' + iD-'^){T + iD"'^) 



a/2 



(17) 



To arrive to Eq. ^ one has to set r = and r' — t, after 
which the square bracket in Eq. ITTl) reduces to {—iDt)'^, 
and integrate dt exp(— iwt). Characteristic times ti in the 



relation Eq. ([TT)) are arranged unevenly as illustrated 
in Fig. [31 The central interval is « t, so that ti are 
located in the close proximity, ti or T2 (see Fig. [3]) either 
to or to t. It is important that in superconducting 
case the arrangement remains the same, and moreover, 
as we will see, Ari and At2 are always <ti 1. This means 
that G{ti — ti+i) can be replaced by Go{ti — ^i+i) times 
the unit matrix. As a result, the matrix structure of 
V drops out. The only Green function that retains the 
matrix structure is G(ti + T2 — t), Fig. [31 However, in 
the component Ln, the anomalous Green function drops 
out, so that 



1*1 



dri dT2 

(TiT2)l-"/2 



G(-ri -T2- t), 



^l+T2<|t 



(18) 

where G{t) is defined by Eq. (fT^ . Eq. dH) immedi- 
ately follows from Eqs. and Eq. dH]). The Green 
function G in Eq. (fTSjl generates the density of states, 
in Eq. One point should be clarified with re- 
gard to the validity of the above result Eq. p^ . We 
used the normal-metal solution Eq. (jl7p . This is jus- 
tified since integrals over ri, T2 in Eq. (|18p come from 
Ti, T2 ~ A~^ exp(— l/a). This also validates the assump- 
tion Ari, At2 <^ 1, which we used to disregard the matrix 
structure of G(ti — i^+i). 



C. Unconventional arrangements of times 

There still remains a question whether or not the ma- 
trix structure of the superconducting Green functions, 
which becomes important near the threshold (w — A) ^ 
A, gives rise to the contributions to A{llj), caused by "un- 
conventional" arrangements of times, ti, {\ti \ ^ A~^), as 



(a) 



(b) 



(c) 







FIG. 4: (Color online) Examples of " unconventional" time do- 
mains in the integrand of Eq. (|16p : (a): Position of the point, 
t\, such that ^ A~^, |i — ti| ^ A~^, does not contribute 
to Ln by virtue of Eq. ((22]); (b); As long as jtil > A"\ 
\t-ti\ > A~\ and \t2\ > A~\ \t-t2\ > A~\ contribution 
of the arrangement of times vanishes upon integration over t\ 
or t2, see Eq. (|26p : (c): For the same reason, "long" (^ A~^) 
intervals in the general "unconventional" arrangement yield 
vanishing contribution to and thus to the absorption at 
the threshold, (a; - A) < A. 
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t t„ t„-l tn-2 ' ' ' tl 



FIG. 5: For exchange interaction with the core hole, "un- 
conventional" arrangement of times, {ti — ti+i) ^ A~^, dom- 
inates the near-threshold, (ctJ — A) ^ A absorption. Odd n 
describes the absorption peak at = A — eo- 

shown in Fig. and Fig. ^p; these arrangements are 
not relevant in the normal-metal case. For example, the 
simplest such "unconventional" arrangement, Fig. 3^, 
manifests itself as an extra combination 

J^dtkG{tk-i - tk)VG{tk - tk+i) (19) 

in the integrand Eq. p^ . Since the arguments of G in 
Eq. are large, one can use the long-time asymptote 

Git) ^Gs{t)( \. "T^'^, (20) 

' ' A|i|»i ""-'y-sgnit) 1 J' 



where the Gs{t) is the At » 1 asymptote of Eq. (fT2|) 

Gs{t)=^oSgn{t) (^|^y^'^e-^l'l+^W4. (31) 

Note however, that the matrix structure in the integrand 
of Eq. mi) is 

(-\-O(J-0(-^")^»' 

Thus, we turn to the next possible arrangement of times 
Fig. [3b; the corresponding combination in Eq. (jl6p com- 
ing from this arrangement reads 

dtk / dtk+iG{tk-i~tk)VG{tk-tk+i)VG{tk+i-tk+2)- 

(23) 

To integrate over t^, we perform multiplication of the 
first three matrices and obtain 



Gs{tk-l—tk] 



-1 1 



1 VGitk-h+i) F{tk-tk+i) 
-l)\F{tk-tk+i) G{tk~tk+i) 



r 



1 1 



= Gs{tk-i-tk)[G{tk-tk+i)+F{tk-tk+i)\ , ^ _^ , 

m 



Then the integration over tk in Eq. (f23]l reduces to 

j dT[G{T)+F{T)\Gs{tu-i-tk+i-T), (25) 



where we introduced a variable r = tk—tk+i- Typical dis- 
tance between the points, \tk-i —tk+i\ and \tk+i —tk+2\, 
is ^ A~^, which suggests that the limits of integration 
can be extended to ±00. Upon this extension we get 



giA(tfc+i-tfc_i) 
\/\tk+l — tk-l\ J- 



dT[GiT) + FiT)]e^^\ (26) 



which is identical zero. The same reasoning rules out^^ 
the more complex "unconventional" arrangements of 
times at the threshold, like the ones shown in Fig. [4);. 
These arrangements, however, become essential in the 
case of exchange interaction with core hole, to which we 
now turn. 



III. EXCHANGE INTERACTION WITH CORE 
HOLE 

Exchange interaction with core hole corresponds to re- 
placement 



where S is a localized spin, and cr is electron spin op- 
erator. To illustrate the dramatic impact which the ex- 
change interaction has on the near-threshold absorption, 
we return to Fig. |4^ and corresponding expression Eq. 
(fT9|) . For potential interaction with core hole, this ex- 
pression was identical zero by virtue of relation Eq. (|22p . 
Recall now that in the stationary problem the diagonal 
part of the exchange interaction, F(r)5^cr^, creates two 
in-gap bound states^^: one below the upper edge by 



(28) 



V(v) ^ Jd{r) (S • a) , 



(27) 



and one above the lower edge by Sq- The reason behind 
this effect is that S^a^ effectively transforms the operator 
V in Eq. (jlOp into the unity matrix. An immediate con- 
sequence of this transformation for our calculation is that 
the contribution Eq. (jl9[) becomes finite. Subsequently, 
the contribution Fig. [Us and all higher-order "uncon- 
ventional" contributions illustrated in Fig. [5^ are also 
finite. Within our formalism, the in-gap bound states 
emerge as poles, ± (A — Eq)], of the Green function 
upon summation^" of infinite series of diagrams. 

In deriving Eq. ([5]) for A{uj) near the threshold, we 
in fact repeat all the steps which would render the sta- 
tionary in-gap states. Namely, we notice that the phase 
^X^fc \tk+i — tk\ of the integrand in Eq. is large, 

which insures that the dominant contribution to Ln(t) 
comes from the domain < ti < t2 ■ ■ ■ < t, see Fig. 
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[5^, when the net phase is At; contributions from the 
domains where tm, are not ordered are suppressed by os- 
cillations of the integrand. Thus we conclude that the 
integral Eq. is dominated by t„i ~ t{m/n). For the 
asymptote Eq. pT|) to be applicable in this domain, the 
condition {tm+i — tm) ^ t/n ^ A^^ must be met. With 
tm ordered, the n-fold integration in Eq. (|16p can be car- 
ried out with the help of the identity 



dx 



a \/[x — a){b — x) 



(29) 



Depending on the parity of n, the remaining integration, 
upon introducing the variables Zi = ti/t, reduces to 



11 1 
dzi j dz2--- j dz(n-l)/2^ 

Zl Z(n-3)/2 



r(^) 



(30) 



for odd n, or to 

1 1 
dzi / dz2 ■ ■ 







for even n. Finally we get 

^nWl, =(-ir 



r(H±i) 



(31) 



iAt 



2 J r(^) 



(32) 



The product, a'^Ln{t), has a sharp maximum at n ~ 
a^At, so that At/n ~ is large, which justifies the 

above assumption {tm+i — tm) ^ A^^. 

The sum over even n, Levenit) = I]euen(o^/2i'o)"in(i), 
leads to the result Eq. Most conveniently it can be 
seen by transforming to the frequency domain, since the 
expansion of Eq. (O in powers of has a form 



This expansion coincides term by term with the sum, 

2p M 



(33) 



L2p{t) exp(-zwi), (34) 



with L2p{t) given by Eq. (|32p. The sum over odd terms 
results in a simple exponent. 



IV. INELASTIC ABSORPTION 



Up to now we neglected the spin- flip part. 



J5{r)[S+a- + S-a+], 



(36) 



of the exchange interaction. As it was mentioned in 
the Introduction, this spin-flip part of interaction be- 
tween electron and core hole creates an effective electron- 
electron scatteringiS.. This explains the possibility of 
inelastic processes with three quasiparticles in the final 
state, as illustrated in Fig. [T|d. The threshold of in- 
elastic process is a; = 3 A. Here we will restrict ourself 
only to the behavior of inelastic absorption away from 
the threshold, (w — 3A) ^ Eq, and follow the calculation 
in Ref.^'*. A great simplification away from threshold is 
that a "golden- rule"- based calculation is sufficient. The 
rate of the process depicted in Fig. [T]d is given by the 
following sum over the quasiparticle states with energies, 
e, e+, and e-. 



£++€_), (37) 



where the first factor is the square of the amplitude, 
which is non-zero since the process involves a spin-fiipi^, 
and the dimensionless spin-flip coupling constant is 



(38) 



Near the threshold, w = 3A, we have e ~ A, e+ « A, 
and e_ « —A. The matrix element near the threshold is 
approximately constant. This simplifies the summation 
in Eq. §7}i to 



W{lu) 



2A2 



(39) 



X deiy{e) jde+i'{e^) / de_i^(e_) (5(a; — e — e+ + e_). 



A — oo 

2A 



Note that in the close vicinity of the threshold, jw — 3A| < 
Eq, in-gap states created by the spinful core hole partici- 
pate in the absorption, as illustrated in Fig. [TJd. Namely, 
a pair of quasiparticles in the final state can consist, e.g., 
of one quasiparticle excited above the gap and empty 
lower in-gap state. 



DISCUSSION 



^odrf(t) f-^j L„(t) cxexp[i(A-eo)i]. (35) 

odd 



This exponent gives rise to the 5-peak, Eq. 
absorption spectrum. 



in the 



Our results Eqs. ([5]) establish the threshold be- 

havior of A[uj) for a general situation when the density 
of states is strongly modified near the Fermi level but 
assumes a constant value away from the Fermi level. A 
notable example is a ID interacting system. The shape of 
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the Fermi-edge singularity in ID interacting electron gas 
in the Luttinger-liquid regime has been studied in^l us- 
ing the bosonization technique. Backscattering plays an 
important role in the exponent of the absorption. When 
backscattering opens a gap, the physics described in the 
present paper comes into play. The case of ID Mott in- 
sulator near half filling makes the behavior of A{lli) even 
richer, since the doping shifts the threshold. A related 
example is the Peierls insulator, when the charge den- 
sity wave and ensuing gap at the Fermi level are due 
to electron-phonon interactions. Note, that in the lat- 
ter case the gap is orders of magnitude larger than in 
superconductor. 

Speaking about conventional setting for Fermi-edge ab- 
sorption in metals, singularity in A{uj) is smeared due to 
the finite lifetime, 7, of the core hole. In our considera- 
tion we assumed that the gap, 2A, exceeds 7. In most 
experiments in metals the smearing of the edge is a frac- 
tion of eV, i.e., much bigger than a typical 2A-value. 
However, the origin of this smearing is not a natural core 
hole lifetime broadening but rather a finite instrumental 
resolution^^. The fact that observed absorption shape is 
a convolution of the singular a Gaussian, which is 

measurement-related, and a Lorentzian, describing nat- 
ural core hole lifetime, allows to separate the two contri- 



butions to the edge smearing. Early attempts^s of such 
separation yielded 7 = 40 meV for 2p core hole. In the 
other experiment^ involving core hole four times shal- 
lower than in Ref. [23l , the natural width was found to be 
four times smaller, 7 = 10 meV. In later experiment—, 
where the full broadening, 29 meV, was very small, anal- 
ysis of the data for the same absorption line as in Ref. [l^ 
revealed even smaller value of the core hole width in sim- 
ple metals, 7 = 4 meV. 

As a final remark, the relevance of the exchange inter- 
action of electron with the core hole was first pointed out 
in Ref. M- 
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